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ABSTRACT 



Aims. We explore the effect of an inhomogeneous mass density field on frequencies and wave profiles of torsional Alfven oscillations 
in solar coronal loops. 

Methods. Dispersion relations for torsional oscillations are derived analytically in limits of weak and strong inhomogeneities. These 
analytical results are verified by numerical solutions, which are valid for a wide range of inhomogeneity strength. 
Results. It is shown that the inhomogeneous mass density field leads to the reduction of a wave frequency of torsional oscillations, in 
comparison to that of estimated from mass density at the loop apex. This frequency reduction results from the decrease of an average 
Alfven speed as far as the inhomogeneous loop is denser at its footpoints. The derived dispersion relations and wave profiles are 
important for potential observations of torsional oscillations which result in periodic variations of spectral line widths. 
Conclusions. Torsional oscillations offer an additional powerful tool for a development of coronal seismology. 
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1. Introduction 

Recent space-based observations revealed a presence of vari- 
ous kinds of magnetohydrodynamic (MHD) waves and oscilla- 
tions in the solar corona. These observations as well as mod- 
eling of MHD waves are important as these waves contribute 
to the coronal heating problem (Roberts 2000) and they may 
consist unique tool of a coronal seismology (Edwin & Roberts 
I1983I Nakariakov & Ofman 2001). Fast kink (Aschwanden et 
al. [19991 Nakariakov et al. [19991 Wang & So lanki 120041 1 and 
sausage (Nakariakov 2003, Pascoe et al. 120071 ) as well as slow 
(de Moortel et al. 2002, Wang et al. 2003) magnetosonic oscil- 
lations were observed to be associated either with or without 
a solar flare. Analytical studies of these oscillations in coronal 
loops were carried on over the last few decades, amongst others, 
by Edwin & Roberts (1982, 1983), Poedts & Boynton (1996), 
Nakariakov (2003), Van Doorsselaere et al. (2004a,b), Ofman 
(2005), Verwichte et al. (2006) and Diaz et al. (2006). 

Coronal loops act as natural wave guides for magnetosonic 
and torsional Alfven waves. The later are purely azimuthal os- 
cillations in cylindrical geometry. In the linear regime, Alfven 
oscillations do not lead to mass density perturbations. As a re- 
sult, contrary to magnetosonic waves, torsional Alfven waves 
can be observed only spectroscopically. While propagating from 
the base of the solar corona along open magnetic field lines, 
these waves may lead to an increase of a spectral line width with 
heigh t (Hassler et al. 119901 Banerjee et al. 119981 Doyle et al. 
119981 ). In closed magnetic field structures, such as coronal loops, 
these waves can be observed indirectly as periodic variations of 
non-thermal broadening of spectral lines (Zaq arashvili 120031 ) . 

Alongside magnetosonic waves, torsional oscillations can be 
used to infer, in the framework of coronal seismology, plasma 
properties inside oscillating loops. These oscillations are an ideal 
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tool of coronal seismology as their phase speed depends alone on 
plasma quantities within the loop, while wave speeds of magne- 
tosonic oscillations are influenced by plasma conditions in the 
ambient medium. Having known mass density within a loop, 
coronal seismology, that is based on torsional oscillations, en- 
ables to estimate a magnetic field strength. Torsional oscilla- 
tions are potentially important in the context of rapid attenua- 
tion of coronal loop kink oscillations (Aschwanden et al. 119991 
Nakariakov et al. 1999). One of a few suggested mechanisms 
of the attenuation is a resonant absorption of fast magnetosonic 
kink waves by azimuthal Alfven waves (Ruderman & Roberts 
2002). This process may lead to a formation of torsional oscil- 
lations in the outer part of a loop. As a result, spotting torsional 
oscillations after the kink mode was attenuated would serve as 
an evidence of this attenuation mechanism. 

A theoretical study of Alfven oscillations in a coronal loop 
was carried on recently by Gruszecki et al. (2007) who con- 
sidered impulsively generated oscillations in two-dimensional 
straight and curved magnetic field topologies. They found that 
lateral leakage of Alfven waves into the ambient corona is neg- 
ligibly small. However, mass density profiles were adopted ho- 
mogeneous within the loop, while the real conditions there are 
much more complex. 

Despite of significant achievements in a development of re- 
alistic models there is still much more effort required to develop 
our knowledge of wave phenomena in coronal loops. A goal of 
this paper is to study the influence of inhomogeneous mass den- 
sity fields on spectrum of torsional oscillations. The paper is or- 
ganized as follows. Analytical solutions for torsional oscillations 
in a longitudinally inhomogeneous coronal loop are presented in 
Sect. 2. The numerical results are showed in Sect. [3] Guidelines 
for potential observations of these oscillations are presented in 
Sect. 4. This paper is concluded by a discussion and a short sum- 
mary of the main results in Sect. [5] 
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2. Analytical model of torsional oscillations 



We consider a coronal loop of its inhomogeneous mass density 
go(z) and length 2L, that is embedded in a uniform magnetic field 
B = Bqz. Small amplitude torsional Alfven waves in a cylindri- 
cal coordinate system (r, <f>, z), in which plasma profiles depend 
on a longitudinal coordinate z only, can be described by the fol- 
lowing linear equations: 



da,/, 
~~dt 



Attqq{z) dz 



(1) 

(2) 



where u$ and are the velocity and magnetic field components 
of Alfven waves. 

These equations can be easily cast into a single wave equa- 
tion 



<9 2 Mrf 
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Fig. 1. Spatial profile of the background mass density, Qo(z), given by 
Eq. l[6} with a 2 = 50. The mass density and length are expressed in 
units of 1CT 12 kg irT 3 and 1 Mm, respectively. 



where Va(z) = Bo/ y/^ngoiz) is the Alfven speed. Assuming that 
Ud, ~ exp(icjt), where to is a wave frequency, we get the equation 



d 2 u 4 , 
~dz T 



+ — M0 = O. 



(4) 



For a trapped solution u$ must satisfy line-tying boundary con- 
ditions which are implemented by setting 



m (z = ±L) = . 



(5) 



Equation with condition (0 consists the well-known Sturm- 
Liuville problem which solution depends on the profile of Va(z)- 
We model the coronal loop by a rarefied plasma at the loop apex 
(at z = 0) and by a compressed plasma at the loop footpoints 
(z = +L). Specifically, we adopt 



Qo(z) = Qoo 1 + a 



2 _ 

L 2 



(6) 



where goo is the mass density at the loop apex and a 2 is a param- 
eter which defines a strength of the inhomogeneity. For a 2 — 
the above mass density profile corresponds to a homogeneous 
loop, while for a larger value of a 2 the medium is more inhomo- 
geneous. FigureQ]illustrates Qo(z) for a 2 = 50. The mass density 
is described by Eq. (O with qqq = 10~ 12 kg rrr 3 and L — 25 Mm. 
Note that plasma is compressed at z — ±L. Substituting Eq. © 
into Eq. ©, we obtain 



d 2 u,p 
"dz 2 



2 1 1 
+ — r- 1 +(T-7 
V 2 A0 I 



0. 



where V A o - Bo/ ^4ngoo- With a use of the notation 

a> L 



y = u 4 



2au> 
V^L Z 



V A o 2a 



(8) 



Eq. (|7]i can be rewritten in the form of Weber (parabolic cylinder) 
equation (Abramowitz & Ste gun [19641 



dx 2 



a\y 



(9) 



Standard solutions to this equation are called Weber (parabolic 
cylinder) functions (Abramowitz & Stegun [T964| ) 
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and yi(x), yo{x) are respectively even and odd solutions to 
Eq. © 



yiix) = 1 +a — 



1\ x 4 
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2. 1 . Two limiting solutions 

Periodic solutions to Eq. (0 can be written analytically in the 
limiting cases: (a) for a large value of a but a moderate value 
of x; (b) for a large x but a moderate a. The first (second) case 
corresponds to a 2 <K 1 (a 2 » 1). 

2.1 .1 . Weakly inhomogeneous plasma 

We consider first the case of a weakly inhomogeneous mass den- 
sity field, i.e. a 2 <K 1. In this case we have 

(7) a < 0, -a » x 2 , p= V^a . (12) 

We adopt the following expansion (Abramowitz & Stegun 
[1964b : 



W(a, x) + iW(a, -x) = 
V2V(a, 0) exp [v r + i(px + tt/4 + v,)] , 
where 



(x/2) 2 | 2(x/2) 4 
' (2p) 2 {2pf 



2/3(x/2) 3 

Vi = o 

2p 



(13) 

(14) 
(15) 
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As a result of relation -a » x 2 we have from Eq. fl3[ 



2.1.2. Strongly inhomogeneous plasma 



W(a,x)= V2>(a,0)exp( — ^-r)cos£ 
W(a, -x) = >/2W(a, 0) exp 



x 

I6p 



f = px + 7T/4 + 



sin^, 
24 p 



The general solution to Eq. (O is 

= c\ W(a, x) + C2 W(a, -x) , 



(16) 
(17) 
(18) 

(19) 



where c\ and ci are constants. 

For a homogeneous loop, i.e. a 2 = 0, we recognize the well 
known solution 



■ c\ cos {kz + n/4) + C2 sin (kz + n/4). 



(20) 



Here wave number k satisfies the following homogeneous dis- 
persion relation: 



OJ 



(21) 



Line-tying boundary conditions of Eq. © lead then to discrete 
values of the wave frequency, viz. 



V 



AO 



2L 1 +a 2 /6 



1,2,3,.. 



(22) 



From this dispersion relation we infer that in a comparison to 
the loop with a homogeneous mass density distribution, poo, the 
weakly inhomogeneous mass density field results in a decrease 
of a wave frequency. This reduction is a consequence of the fact 
that the inhomogeneous loop is denser at its footpoints, so the 
average Alfven speed is decreased. To show this, we compare the 
results for the inhomogeneous loop with the homogeneous loop 
with the same average density, so that both loops contain exactly 
the same mass (Andries et al. 120051 ). We introduce a frequency 
difference 



Aw n = u n - cj„ , 
where 

nn - nn Bq 
w n = — Vao 



2L 2L ^4jrg 

corresponds to the average mass density 

Substituting Eq. d25l l into Eq. d24l >, we obtain 



nn 



(Or, = 



Vao 



2L Vl + a 2 /3 



(23) 



(24) 



(25) 



(26) 



From Eqs. ( f23b and (26) we find that Au> n < 0. Here we infer 
that in comparison to the average mass density case the wave 
frequency is reduced, but as a result of a 2 «: 1 the frequency 
reduction is small. This is in a disagreement with Fermat's law 
and with the results of Murawski et al. ( 120041 ) who showed that 
sound waves experience frequency increase in a case of a space- 
dependent random mass density field. 



We discuss now a strongly inhomogeneous mass density case, 
i.e. a 2 » 1. This case corresponds to x » \a\. In this limit we 
get (Abramowitz & Stegun [T964| l 



W(a, x) = \2k/x(si(a, x) cos(^) - S2(a, x) sin(£)) , 
W(a, -x) = yj2/kx(s\ (a, x) sin(£) - s 2 (a, x) cos(£)) , 

where 

x 2 7T argr(l/2 + id) 

t = — - a In x + - h , 

b 4 4 2 



si(a, x) ~ 1 + 
S2(a, x) ~ 

with 
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r(r + 1/2 + ia)/T(l/2+ia), r = 2,4,.. 



(27) 
(28) 

(29) 

(30) 
(31) 

(32) 
(33) 



The boundary conditions of Eq. (© lead to the discrete fre- 
quency spectrum 



nn V A0 
a L 



(34) 



Here we infer that the strongly inhomogeneous mass den- 
sity field results in a significant decrease of a wave frequency in 
comparison to the case of the loop with the constant density, goo- 
This wave frequency decrease is a consequence of the fact that 
the inhomogeneous loop is denser at its footpoints. Substituting 
Eq. (34) into Eq. d23l we find that Au> n > 0. This wave frequency 
decrease, in a comparison to the case of an average mass density 
is now in an agreement with Fermat's law and with the results of 
Murawski et al. (2004). 

3. Numerical results 

Numerical simulations are performed for Eqs. (Q]i, © with an 
adaptation of CLAWPACK which is a software package de- 
signed to compute numerical solutions to hyperbolic partial dif- 
ferential equations using a wave propagation approach (LeVeque 
2002). The simulation region (-L, L) is covered by an uniform 
grid of 600 numerical cells. We verified by convergence studies 
that this grid does not introduce much numerical diffusion and as 
a result it represents well the simulation region. We set reflect- 
ing boundary conditions at the left and right boundaries of the 
simulation region. 

Figure|2]shows a spatial profile of velocity u</,(z) for a 2 = 50, 
drawn at t — 1000 s (solid line). This spatial profile results from 
the initial Gaussian pulse that was launched at t — in the center 
of the simulation region, at z = 0. It is noteworthy that the sine- 
wave profile of Eq. ( 120) . which is valid for a 2 = (dashed line), 
is distorted by the strong inhomogeneity which takes place for 
the case of a 2 = 50. 

As a consequence of the inhomogeneity wave period is al- 
tered. Figure[3]displays wave period P vs. inhomogeneity param- 
eter a 2 . Diamonds represent the numerical solutions while the 



solid lines correspond to the analytical solution to Eqs. (|22| | (top 
panel) and (|34l l (bottom panel). Wave periods were obtained by 
Fourier analysis of the wave signals that were collected in time 
at the fixed spatial location, z = 0. It is discernible that the nu- 
merical data fits quite well to the analytical curves. A growth of 
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Fig. 2. Numerically evaluated velocity profile u$ at t = 1000 s for a 2 = 
50 (solid line). This profile corresponds to the mode number n = 1 . Note 
that as a result of strong inhomogeneity, u$ departs from the sine-wave 
which corresponds to a 2 = 0. The dashed line corresponds to Eq. (20) 
with c\ = c'2 = 0.5. 



120 




10 20 30 40 50 
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Fig. 3. Wave period P = co/2n vs. a 2 for the mode number n = 1. 
Diamonds correspond to the numerical solutions to Eqs. (TJ, {2^. Solid 
lines are drawn with the use of the analytical solution to Eqs. d22b and 
d34t . The wave period is expressed in seconds. 



wave period P with a- results from wave scattering on centers 
of the inhomogeneity and it can be explained on simple physical 
grounds. In an inhomogeneous field wave frequency co n of the 



torsional oscillations can be estimated from the following for- 
mula: 

nn - 

w n = — V A o, (35) 

where Vao is the averaged Alfven speed that is expressed by Eq. 
(24). Using P = 2nlco a we obtain 

P = \ • (36) 
nB 

As go grows with a, the growth of P with a results in. 



4. Potential observations of torsional oscillations 

Torsional oscillations of a coronal loop may result in periodic 
variations of spectral line non-thermal broadening (expressed by 
a half line width, AAb, hereafter HW) (Zaq arashvili |2003"T > . For a 
homogeneous loop, HW can be expressed as 

uVp&A 

AA B - |sin(w n f)sin(£ n z)| , (37) 

c 

where u is an amplitude of oscillations, A is a wave length of 
the spectral line and c is the light speed. Periodic variations of 
spectral line width depend on a height above the solar surface: 
a strongest variation corresponds to the wave antinode and the 
place of a lack of line width variation corresponds to the nodes 
(loop footpoints). Therefore, time series of spectroscopic ob- 
servations may allow to determine a wave period. Knowing a 
length of the loop, we may estimate the Alfven speed, which in 
turn gives a possibility to infer the magnetic field strength in the 
corona. We estimate the expected value of line width variations 
which result from torsional oscillations. For a typical coronal 
Alfven speed of ~ 800 km/s, an amplitude of linear torsional 
oscillation can be ~ 40 km/s, which consists 5% of the Alfven 
speed. For the "green" coronal line Fe XIV (5303 A) from Eq. 
(37) we obtain 

AAb ~ 0.7 A. (38) 

This value is about twice larger than the original thermal broad- 
ening of Fe XIV line. As a consequence, torsional oscillations 
can be detected in time series of the green coronal line spectra. 

For a weakly inhomogeneous distribution of mass density 
along a loop, Eq. (l22l i enables to estimate the Alfven speed at 
the loop apex with the help of the observed period of HW vari- 
ation and a loop length. For a strongly inhomogeneous density 
profile along a loop, Eq. (1341 shows that a wave period of tor- 
sional oscillations is not just the ratio of the loop length to the 
Alfven speed, but it strongly depends on the rate of inhomo- 
geneity, a 2 . Therefore, an additional effort is required in order 
to apply the method of coronal seismology for torsional oscil- 
lations. A spatial variation of mass density along the loop can 
be estimated by a direct measurement of spectral line intensity 
variation along the loop. Then, the estimated variation can be fit- 
ted to Eq. ©, and hence a value of a 2 can be inferred. Eq. fl34T l 
provides a value of Vao at the loop summit. Another possibility 
is to collect time series of spectroscopic observations at differ- 
ent positions of the loop. A spatial variation of line width along 
the loop may be compared to the theoretical plot of (Fig. 2), 
which enables to estimate a 2 and consequently Alfven speed at 
the loop apex (with a use of Eqs. (l22l i or fl34"li). 
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5. Discussion and summary 

It is commonly believed that Alfven waves are generated in the 
solar interior either by convection (granulation, supergranula- 
tion) or by any other kinds of plasma flow (differential rotation, 
solar global oscillations). Due to their incompressible nature, 
these waves may carry energy from the solar surface to the solar 
corona and therefore they may significantly contribute to coro- 
nal heating and solar wind acceleration. In closed magnetic loops 
the Alfven waves may set up the standing torsional oscillations, 
while in opened magnetic structures these waves may propagate 
up to the solar wind. As a result, observations of Alfven waves 
can be of vital importance to the problems of plasma heating and 
particle acceleration. 

The Alfven waves that propagate along open magnetic field 
lines may lead to a growth of a spectral line width with height 
(Hassler et al. [19901 Banerjee et al. [19981 Doyle et al. [19981 . 
However, at some altitudes the spectral line width reveals a sud- 
den fall off (Harrison et al. 120021 O'Shea et al. 120031 120051 ). This 
phenomenon was recently explained by resonant energy transfer 
into acoustic waves (Zaqarashvili et al. 120061 1. 

On the other hand, the photospheric motions may set up tor- 
sional oscillations in closed magnetic loop systems, which can 
be observed spectroscopically as periodic variations of spectral 
line width (Zaqarashvili 20031. As a result, the observation of 
Alfven waves can be used as an additional powerful tool of coro- 
nal seismology; the observed period and loop mean length en- 
ables to estimate the Alfven speed within a loop, which in turn 
makes it possible to infer a mean magnetic field strength. 

Besides their photospheric origin, torsional Alfven waves 
can be generated in the solar corona in a process of resonant ab- 
sorption of the global oscillations (Ruderman & Roberts 2002 
Goossens et al. 120021 Andries et al. l2005l Terradas et al. l20"06l >. 
These oscillations may excite Alfven waves in the outer inho- 
mogeneous part of a loop, leading to attenuation of global oscil- 
lations and amplification of torsional oscillations. These Alfven 
oscillations can be detected as periodic variations of spectral line 
width. As a consequence, observations of Alfven waves can be 
a key for a determination of a damping mechanism of the loop 
global oscillations. 

Dynamics of torsional Alfven waves in a homogeneous loop 
can be easily solved. However, real coronal loops are longitudi- 
nally inhomogeneous, which leads to alteration of wave dynam- 
ics (Arregui et al.[2005 2007 Van Doorsselaere et al. 2004a,b, 
Donnelly et al. 120061 Dymova & Ruderman 120061 McEwan et 
al. 120061 1. Therefore, the dynamics of Alfven waves in longitudi- 
nally inhomogeneous coronal loops must be understood in order 
to provide analytical basis for potential observations of torsional 
oscillations. 

In this paper we discussed by analytical and numerical 
means evolution of torsional Alfven waves in an inhomogeneous 
mass density field. The analytical efforts resulted in dispersion 
relations which were obtained for a specific choice of an equilib- 
rium mass density profile. These dispersion relations were writ- 
ten explicitly for two limiting cases: (a) weekly inhomogeneous 
and (b) strongly inhomogeneous mass density fields. From these 
dispersion relations we inferred that the inhomogeneity results in 
a wave frequency reduction in comparison to that of estimated at 
the loop summit. This analytical finding is supported by the nu- 
merical data which reveals that frequency reduction takes place 
outside the region of validity of the analytical approach. As a 
result of that we claim that a reduction of wave frequency is 
ubiquitous for the inhomogeneous mass density field we consid- 
ered. This reduction is a consequence of wave scattering on in- 



homogeneity centers and it results from reduction of the average 
Alfven speed within a coronal loop. This frequency reduction 
has important implications as far as wave observations are con- 
cerned. The analytical formulae can be used for estimation of 
coronal plasma parameters and therefore torsional Alfven waves 
consist an additional powerful tool of coronal seismology. 
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